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Abstract. While collisionally excited lines in H II regions allow one to easily probe the chemical composition of the interstellar 
medium in galaxies, the possible presence of important temperature fluctuations casts some doubt on the derived abundances. To 
provide new insights into this question, we have carried out a detailed study of a giant H II region, H 1013, located in the galaxy 
Ml 01, for which many observational data exist and which has been claimed to harbour temperature fluctuations at a level of 
t 2 = 0.03 - 0.06. We have first complemented the already available optical observational datasets with a mid-infrared spectrum 
obtained with the Spitzer Space Telescope. Combined with optical data, this spectrum provides unprecedented information on 
the temperature structure of this giant H II region. A preliminary analysis based on empirical temperature diagnostics suggests 
that temperature fluctuations should be quite weak. However, only a detailed photoionization analysis taking into account the 
geometry of the object and observing apertures can make a correct use of all the observational data. We have performed such a 
study using the Cloudy _3D package based on the photoionization code Cloudy. We have been able to produce photoionization 
models constrained by the observed H/3 surface brightness distribution and by the known properties of the ionizing stellar 
population than can account for most of the line ratios within their uncertainties. Since the observational constraints are both 
strong and numerous, this argues against the presence of significant temperature fluctuations in H 1013. The oxygen abundance 
of our best model is 12 + log O/H = 8.57, as opposed to the values of 8.73 and 8.93 advocated by Esteban et al. (2009) and 
Bresolin (2007), respectively, based on the significant temperature fluctuations they derived. 

However, our model is not able to reproduce the intensities of the oxygen recombination lines observed by Esteban et al. 
(2009), as well as the very low Balmer jump temperature inferred by Bresolin (2007). We have argued that the latter might 
be in error, due to observational difficulties. On the other hand, the discrepancy between model and observation as regards 
the recombination lines cannot be attributed to observational uncertainties and requires an explanation other than temperature 
fluctuations. 
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1. Introduction 

Giant H II regions are among the most effective probes of the 
chemical composition of the present-day interstellar medium in 
galaxies. Methods to derive abundances in these objects have 
been devised many years ago (Aller & Menzel 1945) and are 
straightforward. If the electron temperature can be measured 
using a temperature indicator such as the emission line ratio 
[Om] ,14363/5007 or [Nil] ,15755/6584, the ionic abundances 
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can be obtained directly from the intensities of lines emitted by 
the ions. Elemental abundances are then obtained by correcting 
for unseen ionic stages, as explained e.g. in Osterbrock (1989) 
or Stasinska (2004). 

There are, however, a few issues which have appeared over 
the years for which no broadly accepted solution yet exists and 
which cast some doubts on the derived abundances. One of 
them is the question of the possible presence of important "tem- 
perature fluctuations" in H II regions (which should better be 
called "temperature inhomogeneities"). It has been argued by 



2 



Stasinska et al.: H 1013 



Peimbert (1967) that such fluctuations are likely to exist, since 
the temperature derived from the observed Balmer jump is 
significantly lower than that derived from [O III] /14363/5007. 
Peimbert & Costero (1969) showed that such temperature fluc- 
tuations would bias the derived abundances with respect to 
hydrogen towards values that are systematically too low, un- 
less one adopts their proposed temperature fluctuation scheme 
to account for them. However, up to now, no mechanism has 
been found to produce temperature fluctuations to a level of 
t 2 ~ 0.04 which is a typical value proposed by Peimbert & 
coworkers (Peimbert 1967, Peimbert & Peimbert 2011). It is, 
however, possible that several mechanisms proposed by vari- 
ous authors (e.g. such as those summarized by Torres-Peimbert 
& Peimbert 2003) may add up to produce the observed level 
of temperature fluctuations. Another, perhaps related question, 
is that ionic abundances derived from recombination lines are 
systematically higher than those derived from optical collision- 
ally excited lines (see e.g. Garcia-Rojas & Esteban 2007 for a 
review and a general discussion). These discrepancies could be 
due to the same temperature inhomogeneities that explain the 
difference between the temperatures derived from the Balmer 
jump and the temperatures derived from [O III] /14363/5007. In 
that case, it is believed that abundances obtained from recom- 
bination lines are the correct ones, since recombination lines 
have roughly the same - weak - temperature dependence, while 
optical collisionally excited lines have a very strong temper- 
ature dependence and are thus affected by temperature inho- 
mogeneities (see e.g. Stasinska 2009). Other interpretations of 
these abundance discrepancies are also possible, such as abun- 
dance inhomogeneities (Tsamis & Pequignot 2005, Stasinska 
et al. 2007), density inhomogeneities (Mesa-Delgado et al. 
2012), departure from a Maxwellian energy distribution of the 
electrons (Nicholls et al. 2012), or simply inaccurate theoretical 
rates for the recombination coefficients (Pradhan et al. 201 1). 
Finally, as shown by Stasinska (2005), even in the absence of 
abundance or density inhomogeneities, a further issue is repre- 
sented by the fact that metal-rich H II regions do not appear as 
such when using classical temperature-based abundance deter- 
minations from optical lines, due to an important internal tem- 
perature gradient caused by strong emission of infrared lines in 
the central zones. 

In this paper, we perform a detailed analysis of the extra- 
galactic Hll region H 1013 in the spiral galaxy M 101, in or- 
der to look for some clues to the problem enunciated above. 
H 1013 is a giant Hll region, for which deep spectra have re- 
cently been obtained on very large telescopes (Bresolin 2007 
and Esteban et al. 2009). The oxygen abundances derived by 
these authors from optical collisionally excited lines, with- 
out accounting for temperature fluctuations, are 12 + log(0/H) 
= 8.52 and 8.45, respectively. The abundance of ++ derived 
by Esteban et al. (2009) from recombination lines is larger 
than that derived from [O III] /15007 by a factor of 2.5. From 
this discrepancy, and attributing it to temperature fluctuations, 
these authors find t 2 = 0.037 ± 0.020. They also obtain an es- 
timate of t 2 using a maximum likelihood method from various 
helium line intensities and obtained t 2 = 0.031 ± 0.017. On 
the other hand, from his estimate of the Balmer jump temper- 
ature, Tbj = 5000 ± 800 K, and of the temperature deduced 




Fig. 1. Portion of the Ha continuum- subtracted image of M 101 
presented by Cedres & Cepa (2002) showing the H II region 
H 1013. The apertures corresponding to the optical and infrared 
spectra used in this study are also indicated for reference: the 
long narrow one corresponds to the LRIS data, the short narrow 
one to the HIRES data, and the wide one to the mid-infrared 
observations (note that the slit orientation is arbitrary) . 

from the [O in] ,14363/5007 line ratio, T[0 in] = 7700+250 K, 
Bresolin (2007) obtains a temperature fluctuation factor t 2 = 
0.06+0.02, which implies an oxygen abundance larger by about 
a factor of 4 using the Peimbert & Costero (1969) scheme. 
H 1013 is thus an object where previous studies argued for the 
presence of important temperature inhomogeneities. In order to 
acquire additional constraints, we obtained mid-infrared spec- 
troscopic observations with the Spitzer Space Telescope. Mid- 
infrared lines are crucial in the diagnostic of the physical condi- 
tions of this object, since their emissivities depend very little on 
the electron temperature. In addition, they allow one to resolve 
the above-mentioned degeneracy in abundance determinations. 
In order to make the best use of this new information, we con- 
struct a photoionization model of H 1013 that is as realistic as 
possible and constrained by all the available data, including the 
size and positions of the observing apertures. 

The organization of the paper is as follows. In Section [2] we 
present the observational data used in this study. In Section [3] 
we deduce from them some properties of the ionizing radiation 
field and of the gas density distributions that will be used as 
an input in the photoionization modeling. In Section [4] we de- 
scribe our photoionization modeling procedure. In Section [5j 
we present our results. A brief summary and general conclu- 
sions are presented in Section [6] 

2. Observational dataset 

2.1. Narrow-band filter images 

Cedres & Cepa (2002) used CDD observations in several 
narrow-band filters to compile a catalogue of 338 H II regions 
in the inner parts of M 101 (NGC 5457), also providing infor- 
mation about their fluxes, extinctions, equivalent widths, spa- 
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tial distribution, excitations, radiation hardness, ionization pa- 
rameters and metallicities. H 1013 is identified as the Hll re- 
gion number 299 in their catalogue. 

We use the Ha and Hfi continuum- subtracted images 
(kindly provided by B. Cedres) in our study. These images were 
obtained at the Nordic Optical Telescope with the ALFOSC in- 
strument in direct imaging mode (spatial resolution of 0.189 
arcsec/pix). A description of the reduction process (also in- 
cluding flux calibration and adjacent continuum subtraction) 
can be found in the original paper by Cedres & Cepa. Images 
through additional narrow-band filters (H/5, [O II], [O III], [S II] 
and [S ill]) are also available, however the error bars on the 
integrated line ratios are too large (>50%) to provide useful 
constraints to our models. 

Figure [Upresents a small portion of the original continuum- 
subtracted Ha image, centered on H 1013. The size and loca- 
tion of the apertures used to obtain the optical and infrared 
spectra described in Sects. 12.21 and 12.31 are also indicated for 
reference. As can be noticed, the three apertures cover differ- 
ent portions of the H II region, and none covers the whole re- 
gion. This is an important point that must be taken into account 
when comparing the results of our photoionization model with 
the various spectroscopic observations (or when combining op- 
tical with mid-IR lines). 

The continuum- subtracted Hf3 image is used in Sect. l3.3l to 
obtain the nebular density profile. In Sect. I3.1l we also use the 
total extinction-corrected H/3 flux provided by Cedres & Cepa 
to constrain the properties of the ionizing stellar population of 
H1013. 



H1013 




5650 



5700 



5750 



5800 



X (A) 



Fig. 2. The WR features in the optical spectrum of Bresolin 
(2007). 



2.2. Optical spectra 

Two sets of optical spectra are available. The first one is from 
Bresolin (2007), and was obtained with the Low Resolution 
Imaging Spectrometer (LRIS) spectrograph at the Keck I tele- 
scope. The slit size was 1.5 arcsec x 14 arcsec. The seeing 
during the observations was 0.7 arcsec. Three spectra were 
taken: a blue spectrum: 3300-5600 A (FWHM spectral resolu- 
tion ~5 A); red spectrum: 4960-6670 A (spectral resolution ~ 
4 A); NIR spectrum: 6050-9800 A (spectral resolution ~ 9 A). 
We adopt the dereddened line intensities from Bresolin (2007). 

The second spectrum we consider is that by Esteban et 
al (2009). It was taken with the High Resolution Echelle 
Spectrometer (HIRES) at the Keck I telescope. The spectrum 
covers the 3550-7440 A region with a spectral resolution of 
23,000. The slit width was 1.5 arcsec and the extraction was 
performed over the central 5.76 arcsec. The dereddening pro- 
cedure was slightly different from that adopted by Bresolin 
(2007) but, again, we adopted the dereddened line intensities 
from the original paper. The extraction apertures are indicated 
in Fig. Q] for both spectra (note that their orientation is arbi- 
trary). 

As a quick reference, the reddening-corrected line inten- 
sities with respect to H/3 are reported in the third column of 
Table [3J the label in column 2 indicating whether the data are 
from Bresolin (2007: B07) or Esteban et al. (2009: E09). The 
intensities of interest from the two sets of observations agree 



in general very well within the quoted error bars. A notable 
exception is the case of [On] 372740, whose intensity from 
Esteban et al. (2009) is 40% smaller relative to the one ob- 
tained by Bresolin (2007). Such a large difference is not ex- 
pected. Since the nearby high-order B aimer lines in Esteban 
et al. (2009) deviate from the recombination values by over 
two sigmas, we prefer to use the total [O II] 3727+ intensity 
from Bresolin (2007) to constrain the photoionization models. 
However, we consider it safe to use the [O II] /13726/3729 line 
ratio from Esteban et al. (2009) to obtain an additional density 
diagnostic. Some other lines (e.g. [S II] and [Nil]) have inten- 
sities that differ by more than the assigned error bars. This, a 
priori, can be due to the fact that the observing apertures are 
not the same. 

2.3. Infrared spectrum 

A mid-infrared spectrum of H1013 was obtained with the 
Infrared Spectrograph (IRS, Houck et al. 2004) aboard the 
Spitzer Space Telescope on 26 April 2007 (program ID 30205). 
The Short-High module was used to cover the 9.9-19.6 /urn 
wavelength range with a 1 1 .3 x 4.7 arcsec 2 aperture and a spec- 
tral resolution R ^ 600. We acquired 13 122 s cycles in IRS star- 
ing mode at two nod positions along the slit. The total integra- 
tion time was 53 min. The pipeline data were processed with 

1 Here, and in the remainder of the text, 3727+ stands for 3726 + 
3729. 
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Fig. 3. Spitzer Space Telescope spectrum of H 1013. 

Table 1. Spitzer spectrum line measurements (Slit: 11.3 x 4.7 
arcsec 2 ). 



Line 


AQmi) 


Flux (xlO 14 erg cm 2 s l ) 


I(Hl7-6= 1) 


[Siv] 


10.5 


1.63+0.25 


4.4+1.6 


Hi 7-6 


12.4 


0.37+0.12 


1.0+0.3 


[Ne II] 


12.8 


15.2+1.4 


41+14 


[Ne in] 


15.6 


4.95+0.20 


13+4 


[S m] 


18.7 


12.2+0.4 


33+11 



IRSCLEAN (version 1.9) to create bad pixels masks, and with 
the SPICE GUI software (version 2.0.2) for spectral extraction. 
Since the spectra obtained from the two nod positions agree 
within the statistical noise, we combined them to yield the fi- 
nal, flux-calibrated spectrum of H 1013, shown in Fig. 3. The 
fluxes of the emission lines present in the spectrum were mea- 
sured with the splot task within iRAlU, and are summarized in 
Table 1, both in terms of absolute flux (in erg cm" 2 s" 1 ) and rel- 
ative to the H 1 7-6 line. In principle, the absolute calibration of 
the mid-IR data is correct within 30%, but we prefer to use the 
mid-IR and optical hydrogen lines to rescale the mid-IR spec- 
trum (see later). This ensures that line flux ratios are correct 
across the whole wavelength range. 



3. Derivation of stellar and nebular properties 

3. 1. The ionizing stellar population 

The first thing to note is that the spectrum of H 1013 contains 
Wolf-Rayet (WR) features (see Fig.0. This means that, if the 
ionizing cluster can be approximated by an instantaneous star- 
burst of large enough mass, its age is around 2-5 Myr. 

To obtain the total number of ionizing photons, we as- 
sume in the first place that all the stellar ionizing photons 
are absorbed by the nebular gas. We adopt a distance to 
M 101 of 7 Mpc (a compromise between the value of 7.5 Mpc 
from Sandage & Tammann 1976 and that of 6.85 Mpc from 
Freedman et al. 2001). From a total extinction-corrected flux 



2 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 4. The time dependence of several characteristics of 
an instantaneous starburst of solar metallicity obtained with 
Starburst99 (see text for details). The observed values or their 
upper limits in H 1013 are shown respectively by green and 
magenta horizontal lines. 



in U/3 of 3.6 i 1.2 x 10" 13 erg cm" 2 s" 1 (Cedres & Cepa 
2002), we obtain a total number of hydrogen ionizing photons 
Qn = 4.42 x 10 51 s _1 . Adopting a Kroupa (2001) stellar initial 
mass function and using the Starburst99 stellar synthesis code 
(Leitherer et al. 1999) to relate stellar masses, ages and ioniz- 
ing photon output, we find that, for an age of 2.5 Myr (appro- 
priate for our object, as will be shown below), this corresponds 
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to a total of ~ 1.3 x 10 5 M© of stars with initial masses be- 
tween 0.1 and 100 M©. In fact, dust must be present inside the 
H II region, since iron is heavily depleted (Esteban et al. 2009), 
and therefore absorbs part of the ionizing photons. In our fi- 
nal models (see below), dust absorbs about 1.5 x 10 51 ph s" 1 , 
therefore the total initial stellar mass of the ionizing cluster is 
about ~ 1.7 X 10 5 M . According to Cervino et al. (2003), the 
minimum initial cluster mass necessary to derive trustworthy 
results from photoionization models at an age of about 2.5 Myr 
is at least 10 5 M for the ions we are interested in, if using clas- 
sical synthesis models that assume full sampling of the initial 
mass function. Those authors however recommend a value 30 
times larger for really safe results. Clearly, our object does not 
comply with such conditions. In the absence of direct informa- 
tion on the ionizing stars (such as is sometimes available, see 
Jamet et al. 2004), a Monte-Carlo approach would in principle 
be needed. However, since the hardness of the ionizing radia- 
tion field is directly linked with the observed WR features, we 
may use the stellar population synthesis code Starburst99 to 
estimate the spectral energy distribution that produces the WR 
features observed in H 1013. 

Figure |4] shows the time dependence of several features 
obtained by running Starburst99 with the expanding non-LTE 
stellar atmosphere models implemented by Smith et al. (2002), 
using the Geneva tracks and high mass-loss rates. The run of 
Starburst99 was made for solar metallicity, which is the most 
appropriate for our object. From top to bottom, Fig. 21 shows 
the values of 

1. EW(Ha), the equivalent width of the Ha emission for an 
instantaneous burst, 

2. the ratio QwJQw of helium to hydrogen ionizing photons 
predicted by the model, 

3. the intensity of the stellar Heii/14686 feature (2.5 x 10" 15 
erg cm" 2 s" 1 ), 

4. the intensity of the stellar N in AA6AQ + C in AA65Q feature 
with respect to He II /14686, 

5. the intensity of the stellar C ill A5696 feature with respect 
to Hen ,14686, 

6. and the intensity of the stellar C IV /15808 feature with re- 
spect to He II /14686 . 

The observed values are represented with the green hori- 
zontal lines, the error bars with red dotted lines, and the up- 
per limits (for C in /15696 and C IV /15808) are represented with 
magenta lines. The QnJQn ratio is actually not observed, but 
this value determines the observed value of the emission line 
ratio [O III] A5001/[O II] A3121 in the nebula. From Fig. H we 
see that the observed equivalent width of the Ha line indicates 
an age of 2-3 Myr, the intensity of the stellar He II /14686 fea- 
ture is compatible with ages of about 2.5, 3.2 and 4.2 Myr, 
while the N in /14640 + C in /14650 feature indicates an age of 
2.5-3 Myr, or larger than 5 Myr. On the other hand, the absence 
of observed stellar carbon features indicates an age smaller than 
3 Myr. In fact, in the framework of an instantaneous starburst, 
the observational constraints on the stellar population are al- 
most incompatible, but taking into account uncertainties, they 
point to an age of 2.5 Myr. However, if we consider that, in 
practice, a starburst is never instantaneous but occurs during a 
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Fig. 5. Plasma diagnostics using PyNeb. The dashed or dot- 
ted values indicate the density-temperature loci corresponding 
to the observed reddening-corrected line ratios (also aperture- 
corrected when combining Spitzer and optical data, as ex- 
plained in Sect. I?]). The coloured bands delineate the one sigma 
error zones in the corresponding ratios. Diagnostics involv- 
ing oxygen ions are in green: [Oll]na 3727+/7325+; [Oll]n: 
3729/3726; [Om]na: 5007/4363. Diagnostics involving sulfur 
ions are in black: [Sn]n: 6716/6731; [Sn]na: 6720+/4072+; 
[Sm]fn: 18.7//m/9069; [Sm]na: 9069/6312; the [Nn]na 
6584/5755 diagnostic is in blue; the [Ne m]fn 15.6//m/3869 di- 
agnostic is in pink; finally the [Fe III] 5272/4987 diagnostic is 
in blue. (The characters n, a, f in the denomination of the line 
ratios stand for 'nebular', 'auroral' and 'fine- structure' , respec- 
tively.) 

time interval of a few 10 5 yr, the range of possible solutions 
is larger. We approximate such a situation by considering two 
instantaneous starbursts of different intensities separated by a 
small time interval. After trials and errors, we find that combin- 
ing one burst giving <2h = 4.6xl0 51 ph s" 1 at an age of 2.5x1 6 
yr and another one giving <2h = 1 x 10 50 ph s" 1 at an age of 
3 x 10 6 yr correctly reproduces the observational constraints on 
the stellar populations and allows us to correctly reproduce the 
observed emission line ratios of the nebula (see below). 

3.2. Simple plasma diagnostics 

Before proceeding to the description of the photoionization 
modeling, it is useful to consider the temperature-density 
diagnostic diagram for H1013. This diagram, presented in 
Fig. is obtained by considering all the observed temper- 
ature and/or density sensitive intensity ratios of lines of the 
same ion that have been observed in H1013. We use the 
extinction-corrected intensity measurements and associated 
uncertainties from Bresolin (2007), except in the case of the 



3 This combination of ages does not have to correspond to reality 
since, as explained above, we are in the regime where statistical fluc- 
tuations matter. The important thing is that the spectral energy distri- 
bution of the ionizing radiation field is correct. 
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[O II] /13726/3729 ratio, which was not available in those obser- 
vations and is taken, together with the [Fe III] A 5272/4987 line 
ratio, from Esteban et al. (2009). Note that line ratios involv- 
ing one optical and one infrared line are obtained by forcing 
the H 7-6 / UJ3 ratio to match the case B recombination value 
of 0.0098 for a temperature of 10000 K and a density of 100 
cm" 3 (Storey & Hummer 1995), and the uncertainties are dom- 
inated by the uncertainty in the flux of the H 7-6 line. The dia- 
gram has been constructed with the software PyNeb (Luridiana 
et al. 2012) using the same atomic data as in the photoioniza- 
tion Cloudy (Ferland et al. 1998) in the version adopted for 
the present study. This diagram is of course only indicative, 
since, for each line ratio, it assumes that both density and tem- 
perature are constant. Note that the effects of recombinations 
are not accounted for in such a diagram, while they are taken 
into account in Cloudy. In the present case, this could affect 
the [O II] /13727+/7325+ ratio. Another approximation is that, 
when matching the optical and infrared spectra using the hy- 
drogen lines, one ignores the ionization structure that affects 
the [Ne III] and [S III] lines, which are observed both in the op- 
tical and in the infrared ranges, but through different apertures. 

In addition to the diagnostics obtained by Bresolin (2007) 
and Esteban et al. (2009), [Fein] /15272/4987 provides a use- 
ful diagnostic, as it restricts the density in the corresponding 
emitting zone to a value roughly below 160 cm" 3 . The [Sll] 
A67 16/6731 diagnostic indicates an upper limit of the density 
of 50 cm" 3 in the [S II] emitting zone. With the Spitzer Space 
Telescope observations, we have two additional temperature di- 
agnostics, that could provide some clues on the possible pres- 
ence of large temperature inhomogeneities. If we adopt the 
Peimbert (1967) temperature fluctuation scheme, the values of 
T = 5500 K and t 2 = 0.06 inferred by Bresolin (2007) would 
imply a temperature derived from the [Nelll] /U5.6yum/3869 
ratio of 6150 K, and a temperature derived from the [S III] 
^18.7yum/9069 ratio of 5500 K. A value of t 2 = 0.03 - 0.04, 
as derived by Esteban et al. (2009) would imply a [Nelll] 
^15.6yum/3869 temperature of 6400 - 6300 K, and a [S in] 
^18.7yum/9069 temperature of 6100 - 5900 K. Such values are 
well below the ones obtained from the Spitzer measurements, 
which are roughly 7000 ± 200 K and 9200+*™ K, as inferred 
respectively from Fig. [5] 

While the diagram presented in Fig. [5] is useful to visual- 
ize the various temperature and density diagnostics, it includes 
several approximations, as mentioned above. The only way to 
account for all the observational constraints in the most accu- 
rate way is to build a tailored photoionization model that is 
able to reproduce each of the observational constraints within 
the uncertainties. 



3.3. Global density structure 

The first requirement for a successful photoionization model is 
to be able to reproduce the observed HJ3 surface brightness dis- 
tribution. We use the UJ3 (extinction corrected, continuum sub- 
tracted, and flux calibrated) image from Cedres & Cepa (2002) 
and assume a uniform electron temperature in order to obtain 
the surface brightness distribution in two perpendicular direc- 
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Fig. 6. Upper panel: Density law needed to fit the Hfi surface 
brightness distribution assuming a spherical geometry. Lower 
panel: The observed UJ3 surface brightness distribution in two 
perpendicular directions (black continuous and dotted lines) 
and the computed one assuming a constant electron tempera- 
ture. 



tions. The two distributions are quite similar, indicating that the 
nebula can be assumed to be spherically symmetric. We use this 
HJ3 surface brightness distribution to obtain the density law. We 
assume spherical symmetry and find that a density law of the 
type: 



A^ H =A 1 exp(-(r/A 2 ) A3 ) 



(1) 



is convenient to reproduce the observed H/3 surface brightness 
distribution assuming a constant electron temperature. We find 
a good fit with Ai=13 cm" 3 , A2=107 pc, A3=1.6, as shown in 
Fig. [6] Note that in this case we do not consider any filling fac- 
tor. In the models, the value of A\ is adjusted to reproduce the 
density and ionization diagnostics, and a filling factor smaller 
than unity is then needed to reproduce the observed nebular 
size. In the process of our model fitting, we check that the com- 
puted radial temperature variation has not drastically modified 
the HJ3 surface brightness distribution. 

4. Photoionization modelling procedure 

The photoionization modeling is performed using the code 
Cloudy, version c08.01 (Ferland et al. 1998), within the py- 
Cloud)@ environment. With pyCloudy, we can easily obtain 
the line intensities in specific apertures. The nebula is assumed 
to be spherically symmetric, with the density distribution de- 
scribed in Sect. 13.31 The stellar ionizing radiation field is ob- 
tained as described in Sect. 13.11 The computed line intensities 
are compared with the observed ones, each in its corresponding 
aperture. The effects of seeing for the optical data and of minor 
irregularities in the nebular density distribution are taken into 



4 pyCloudy is the python version of Cloudy _3D (Morisset, 2006) 
available at https://sites.google.com/site/pycloudy/ 
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Table 2. Line ratios used to constrain the model and the asso- 
ciated tolerances. 



AO/O 


line ratio 




density indicators 


0.15 


[OII]3726 / [OII]3729 


0.40 


[Fe III]4702 / [Fe III]4659 


0.40 


[FeIII]4659 / [FeIII]5271 


15 




0.15 


[C1IIIJ5538 /[C1IIIJ5518 




temperature indicators 


0.30 


[NeIII]3869 / [Ne III] 15.5 


0.12 


[OIII]4363 / [OIII]5007 


0.30 


[S III]63 12 / [S III]9069 


0.40 


[S III]9069 / [S III] 18.7 


0.10 


[NII]5755 / [NII]6584 


0.10 


[OII]7325+ / [OII]3727+ 


0.10 


[S II]4070+/[S II]6716+ 




ionization structure 


0.30 


[S III]9069 /[S II]6716+ 


0.40 


[S IV] 10.5 /[S III] 18.7 


0.15 


[O III] 5007 / [Oil] 3727+ 


0.40 


[Ne III] 15.5 / [Nell] 12.8 




abundances 


0.05 


Hel 5876 / Ra 


0.06 


[NII]6584/Hor 


0.05 


[OIII]5007 /HjS 


0.06 


[S II]6716+/Hor 


0.08 


[Ne III] 3869 /HjS 


0.08 


[C1III]5525+ IHJ3 


0.08 


[ArIII]7135 /Ha 


0.15 


[Fe III]4659 / Ha 



account by convolving with a gaussian (adopting a width of 1 
arcsec for the optical lines and 0.5 arcsec for the infrared lines). 

The initial chemical composition of the gas is that derived 
from collisionally excited lines (and taking t 2 = 0) by Bresolin 
(2007) for O, N, and Ne and by Esteban et al. (2009) for He, 
S, CI, Ar, Fe, which were not given by Bresolin. The reason 
why we use the O, Ne, Ne abundances by Bresolin is due to 
the potential problem in the [O II] /13727 intensity of Esteban 
et al. (2009) as explained in Sect. 12.21 The carbon abundance 
is taken equal to the value derived by Bresolin (2007) from the 
C II /14267 recombination line, since there is no other direct 
information on the carbon abundance. 

In the classical logarithmic units taking 12 for hydrogen, 
the initial composition of the gas is thus: H= 12, He= 10.87, 
C= 8.66, N= 7.57, 0= 8.52, Ne= 7.41, S= 6.91, Cl= 5.50, 
Ar= 6.35, Fe= 5.74. 

Since iron is observed to be strongly depleted in this object, 
we consider the presence of dust in the model, using the "grains 
ism" option of Cloudy (but see below for the radial distribution 
of the grains). In the present object, the effect of grains is essen- 
tially to absorb part of the ionizing photons (about 30%) and to 
heat the gas, especially in the central zone. 

A satisfactory model must reproduce the observed YL/3 sur- 
face brightness distribution law (which in our case it does by 
construction), the total reddening-corrected HJ3 flux and the ob- 
served nebular angular size, as well as each of the reddening - 
corrected observed line ratios in each slit, within the observa- 
tional uncertainties. We emphasize that it is not sufficient to 
perform a x 2 fitting of the sum of the deviation of line intensi- 
ties to the observed value: each line carries its own information, 
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and it is important that all observational constraints be properly 
accounted for. 

Absolute calibration of spectroscopic observations is diffi- 
cult. We intercalibrate the Spitzer and optical data by forcing 
the measured value of the H 7-6 / HJ3 ratio to the one predicted 
by our photoionization models in the corresponding slits. The 
value of /(IR), representing the factor by which the measured 
IR fluxes have to be multiplied in order for the H 7-6 / H/3 ra- 
tio to be in agreement with the model, turns out to be -1.07. 
To judge a model, we adopt the same policy as in Morisset & 
Georgiev (2009) and Stasinska et al. (2010). For each observ- 
able O we fix a tolerance r(O), given by 

r(O) = log(l + AO/O), (2) 

where AO/O is the maximum "acceptable" deviation. The val- 
ues adopted for r(O) take into account the observational un- 
certainties and are chosen according to our experience with 
model fitting. For example, while the observational uncertainty 
in [Oin]/i5007/[Oll]/i3727 is small, we take a tolerance of 
1.15 to reflect the fact that this line ratio is very sensitive to the 
density distribution and/or the slit position in our object. 

It is convenient to divide the line ratios to be fitted in dif- 
ferent categories: 

- Ratios of lines measuring the density. 

- Ratios of lines measuring the electron temperature. 

- Ratios of hydrogen lines or of helium lines, which test the 
reddening correction, among others. 

- Ratios of two different ions of the same element, testing the 
ionization structure. 

- Ratios of lines used to determine the chemical composition, 
once the ionization structure is correct. 

The model is compared to the observations by looking at 
the values of 

k(0) = (logO mod - logO obs )/r(0), (3) 

which can easily be plotted in the same diagram for all the line 
ratios. The list of line ratios and associated values of AO/O is 
given in Tabled 

During the fitting procedure, we vary the values of Ai, the 
parameter defining the density (see Sect. 13. 3K the filling factor, 
the total number of ionizing photons (since the proportion of 
those photons that are absorbed by dust is not known a priori), 
and the elemental abundances. We also had to modify the spa- 
tial distribution of the internal dust. Indeed, when considering 
a uniform dust-to-gas ratio in the entire nebula, the inner zones 
became too hot, and the [O III] /14363/5007 ratio could not be 
fitted to gether with the temperature diagnostics of the low ex- 
citation zones. We adopted the simplest model for the dust- to- 
gas ratio, which is a step function. The parameters of the step 
function were obtained by trial and error to reproduce the ob- 
served temperature constraints. Note that such a dust distribu- 
tion which may seem ad hoc is actually fully justified by the 
effects of radiation pressure and the possible melting of grains 
close to the star. That this is precisely what occurs near the 
Trapezium stars in the Orion nebula was shown already long 
ago (Baade & Minkowski 1937). 




Fig. 7. Comparison of model M 1 with the observations. Blue diamonds represent data from Bresolin (2007), green triangles 
data from Esteban et al. (2009). Red symbols represent ratios involving an infrared line. The dashed and dotted lines represent 1 
sigma and 2 sigmas deviations, respectively. 




Fig. 8. Comparison of the composite model, M 2, with the observations. Same layout as in Fig.[7J 



5. Results 

5.1. The simplest model 

In order to constrain the model-fitting we took the observations 
by Bresolin (2007), using the data of Esteban et al. (2009) as a 
complement. The most satisfactory model we produced, within 
the explained procedure, is one having A\ = 200 cm -3 , a fill- 
ing factor of 0.005, and a total number of ionizing photons 
2(H) = 5.9xl0 51 s" 1 . 

A graphical representation of how the line-ratio constraints 
match the data is shown in Fig. [TJ the values of k(0) are plot- 
ted for all the constraints O. Diamonds correspond to inten- 
sities from the Bresolin (2007) spectrum, triangles to intensi- 
ties from the Esteban et al. (2009) optical spectrum, each one 
within its extraction aperture. It can be seen that each of the val- 
ues of k(0) is found to lie between -2 and +2 (except for the 
[S IV] A10.5jnm/[S in] ,11 8.7/rni ratio, and the [O in] ^5007/HyS 
ratio from Bresolin 2007). A more conventional comparison of 
our best model, Ml, with the observed line ratios is shown in 
Table [3j It can be seen that the k(0) values for Bresolin (2007) 
and Esteban et al. (2009) are slightly different, which means 



that either the density distribution in the model is too simplis- 
tic (i.e. that some slight density inhomogeneities are present) 
or that the observational error bars are slightly underestimated. 
Probably both reasons are relevant, but in any case the differ- 
ences are small (except for the case of the [O II] /13727 line 
mentioned above). 

The elemental abundances in this model are: H= 12, He= 
10.97, C= 8.66, N= 7.79, 0= 8.42, Ne= 7.80, S= 6.98, 
Cl= 5.18, Ar= 6.20, Fe= 5.77. They differ from the initial 
ones, mainly due to the fact that the classical ionization correc- 
tion factors used to derive abundances by Bresolin (2007) or 
Esteban et al. (2009) are not in agreement with our photoion- 
ization model. Such an explanation, obviously, does not hold 
for oxygen, because both + and ++ lines are measured. Yet, 
there is a 0.1 dex difference between the model value and the 
value derived by Bresolin (2007). The reason is that this model 
is not perfect, because it cannot reproduce at the same time all 
the line ratios involving [O III] A5001. 

In Fig. [9] we show the radial distribution of the electron tem- 
perature, the dust-to-gas ratio, the hydrogen density, and the 
relative abundances of + and ++ in model Ml (red curves). 
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Table 3. Observed line ratios in H 1013 compared to the results of models Ml and M2. 



line intensities a 


references 


observed 


mod. Ml 


mod. M2 


[0 II] 3726 


E09 


57.00 ± 2.00 


63.02 


91.17 


[0 II] 3729 


E09 


78.00 ± 3.00 


81.51 


114.86 


[Oil] 3727+ 


B07 


221.00+ 11.00 


155.71 


219.82 


[On] 3727+ 


E09 


135.00 + 5.00 


145.16 


206.9 


H13 3734 


E09 


1.50 + 0.10 


2.29 


2.28 


H12 3750 


B07 


3.00 ± 0.20 


2.92 


2.9 


H12 3750 


E09 


2.00 ± 0.20 


2.91 


2.89 


Hll 3771 


B07 


3.30 ± 0.20 


3.81 


3.79 


Hll 3771 


E09 


3.10 + 0.20 


3.80 


3.77 


H10 3798 


B07 


4.40 ± 0.20 


5.11 


5.08 


H10 3798 


E09 


4.20 ± 0.20 


5.10 


5.07 


H9 3835 


B07 


7.30 + 0.30 


7.09 


7.05 


H9 3835 


E09 


6.40 ± 0.30 


7.07 


7.03 


[Nera] 3869 


B07 


3.10 + 0.20 


2.94 


3 


[Nem] 3869 


E09 


2.80 ± 0.20 


3.24 


3.17 


H7 3970 


B07 


16.40 ± 0.70 


15.57 


15.52 


H7 3970 


E09 


15.80 ± 0.60 


15.55 


15.5 


He 1 4026 


B07 


1.54 + 0.08 


1.91 


1.98 


He I 4026 


E09 


1.70 + 0.10 


2.04 


2.12 


[S II] 4072+ 


B07 


1.27 + 0.07 


1.58 


1.58 


[S II] 4072+ 


E09 


1.06 + 0.18 


1.39 


1.39 


H6 


B07 


26.20+ 1.10 


25.48 


25.43 


US 


E09 


25.30 ± 0.90 


25.47 


25.41 


Cn 4267 


B07 


0.25 ± 0.03 


0.30 


0.28 


Cll 4267 


E09 


0.28 ± 0.08 


0.33 


0.31 


Hr 


B07 


46.40+ 1.90 


46.37 


46.32 


Hr 


E09 


44.00+ 1.00 


46.35 


46.3 


[Om] 4363 


B07 


0.24 ± 0.03 


0.20 


0.22 


[Om] 4363 


E09 


0.19 + 0.07 


0.22 


0.23 


He I 4388 


B07 


0.40 ± 0.05 


0.51 


0.53 


He I 4388 


E09 


0.52 ± 0.09 


0.55 


0.57 


He I 4471 


B07 


3.90 ± 0.20 


4.10 


4.26 


On 4651 b 


E09 


0.62 ± 0.27 


0.12 


0.15 


[Fe in] 4659 


B07 


0.40 ± 0.05 


0.34 


0.41 


[Fern] 4659 


E09 


0.25 ± 0.04 


0.32 


0.4 


[Fern] 4702 


E09 


0.14 + 0.03 


0.10 


0.12 


[Fem] 4987 


B07 


0.35 ± 0.04 


0.01 


0.01 


[Fe III] 4987 


E09 


0.32 ± 0.04 


0.01 


0.01 


[Om] 5007 


B07 


102.00 ± 4.00 


84.46 


100.15 


[Ora] 5007 


E09 


97.00 ± 3.00 


93.16 


107.68 


[Fem] 5271 


B07 


0.15 + 0.02 


0.20 


0.24 


[Fem] 5271 


E09 


0.18 + 0.04 


0.19 


0.24 


[CI m] 5518 


B07 


0.46 ± 0.03 


0.41 


0.44 


[CI m] 5518 


E09 


0.39 ± 0.05 


0.42 


0.45 


[CI m] 5538 


B07 


0.29 ± 0.02 


0.29 


0.31 


[CI m] 5538 


E09 


0.25 ± 0.04 


0.30 


0.32 


[Nil] 5755 


B07 


0.59 ± 0.03 


0.51 


0.49 


[Nil] 5755 


E09 


0.43 ± 0.05 


0.44 


0.43 


He I 5876 


B07 


12.10 + 0.60 


11.74 


12.21 


He I 5876 


E09 


11.30 + 0.50 


12.52 


13.07 


[Oi] 6300 


B07 


1.21 + 0.07 


0.16 


0.35 


[0 1] 6300 


E09 


1.44 ± 0.09 


0.12 


0.27 


[S in] 6312 


B07 


0.77 ± 0.05 


1.07 


0.74 


[S in] 6312 


E09 


0.72 ± 0.06 


1.08 


0.74 


Ha 


B07 


293.00 ± 16.00 


294.77 


295.06 


Ha 


E09 


282.00 ± 16.00 


294.71 


295.04 


[Nil] 6584 


B07 


71.00 + 4.00 


73.88 


71.84 


[Nil] 6584 


E09 


59.00 ± 3.00 


64.77 


64.96 


He I 6678 


B07 


3.40 ± 0.20 


3.34 


3.47 


He I 6678 


E09 


3.10 + 0.20 


3.56 


3.72 


[S II] 6716 


B07 


17.80+1.00 


17.25 


16.9 


[S II] 6716 


E09 


14.90 ± 0.90 


15.05 


14.82 


[S II] 6731 


B07 


12.00 ± 0.70 


12.99 


12.92 


[S II] 6731 


E09 


11.00 + 0.70 


11.45 


11.46 


[Arm] 7135 


B07 


7.10 + 0.40 


7.21 


7.05 


[Arm] 7135 


E09 


6.80 ± 0.50 


7.71 


7.56 


[On] 7325+ 


B07 


2.40 ± 0.20 


1.85 


2.65 


[S m] 9069 


B07 


24.00 ± 5.00 


33.10 


23.42 


[SiV] 10.51yum 


IR 


4.40+ 1.94 


2.81 


6.2 


[Neil] 12.81/mi 


IR 


41.00+ 17.01 


30.92 


32.84 


[Nem] 15.55/mi 


IR 


13.00 ± 4.86 


7.45 


9.93 


[S in] 18.67yum 


IR 


33.00+ 13.36 


56.81 


41.02 


log L(Ho') [erg s _1 ] 




39.78 


39.78 


39.78 



a In units of H/? = 100 for optical lines and of H7-6 12.4/mi= 1 for infrared lines 
b Represents the sum of the multiplet 1 of O II 



We see that the radial gradient of the electron temperature is ^ 7700 K in the outer zone. This is what we would have ex 
mild, with a temperature of ^ 7200 K in the central zone and pected for an H II region of roughly solar metallicity. The im 
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Fig. 9. Radial distribution of the electron temperature, the dust- 
to-gas ratio, the hydrogen density, and the relative abundances 
of + (dashed curve) and ++ (continuous curve). The thick 
red curves correspond to Model Ml (or the high-density com- 
ponent of model M2) while the thin blue curves correspond to 
the low density component of model 2. 



portant point is that this temperature distribution is fully com- 
patible with the line intensities observed both in the optical and 
in the infrared. The small increase of the electron temperature 
at 2 x 10 20 cm is due to photoelectric heating of the gas by dust 
grains. We have checked that the amount of dust present in the 
model is compatible with the extinction C(Hj3) ^0.3 estimated 
for this object from the B aimer lines. 



We note that the intensity of the [O I] /16300 line is smaller 
by a factor -10 compared to the observations. This does not 
represent a real worry, since [O I] is essentially produced in 
the ionization front (i.e. not in the ionized regions as all the 
remaining lines), and we did not attempt to model this zone in 
detail. 

The discussion of the recombination lines is postponed to 
SectO 

5.2. A composite model 

One problem that still needs examination is the density struc- 
ture. Our best model, Ml, requires a volume filling factor 
of 0.005. This is not uncommon for giant Hll regions (eg. 
Luridiana et al 1999; Stasinska & Schaerer 1999). Of course, 
using a filling factor is an approximation. There has to be some 
diffuse gas between the filaments. We have to check whether 
such gas affects the intensities of the observed emission lines. 
We have thus computed additional low-density models with the 
same abundances at the "high" density, main model. The com- 
putations are performed until the outer radius of the nebula is 
reached (so, for low densities, these additional models are den- 
sity bounded). A composite model is then obtained by adding 
the fluxes of the main, "high" density model, with one of those 
additional models. This is, of course, a simplification of a more 
realistic model where the density distribution would be fila- 
mentary (as in the model of Jamet et al. 2006), but it is suffi- 
cient for our purposes. The results of our "best fit" composite 
model, M2, are shown in Fig. [8] and Table [3] The high- and 
low-density components have A\ = 235 and A\ = 5 cm -3 , and 
filling factors of 0.003 and 0.5, respectively. 

Characteristics of the corresponding low-density model 
are plotted in Fig. |9j with blue curves. From Fig. [8] we see 
that our best composite model fits the data even better than 
model Ml. In particular, the [S IV] 10.5 / [S III] 18.7 and the 
[O III] /15OO7/H0 ratio from Bresolin (2007) are now well fit- 
ted. This is not too surprising, since we have allowed more 
parameters to vary and the resulting model should be closer 
to reality. Since, for the composite model, we have adopted 
the same dust-to-gas ratio distribution as for model Ml, the 
jump in electron temperature at the place where grains appear 
is much higher than in model Ml, because, as a consequence of 
the higher local ionization parameter due to the lower gas den- 
sity, the effect of heating by grains is enhanced (Stasinska & 
Szczerba 2001). Note that, in the central zone, devoid of grains, 
the temperature of the low-density model is lower than in the 
case of model Ml. This is because oxygen is entirely in the 
form of ++ , the most efficient coolant. In spite of the signif- 
icant differences in temperatures with model Ml, the overall 
"temperature fluctuation" t 2 (as defined by Peimbert 1967 and 
computed by us integrating the temperature weighted by the 
product n e riH over the emitting volume) is very small in this 
composite medel (less than 0.001). We have to note that, al- 
though other composite models can probably be found to fit all 
the observations equally well, they will not differ significantly 
from M2 because of the very strong constraints imposed on 
each zone of the nebula by the observations. 
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The chemical composition of the composite model is : H= 
12, He= 1 1.0, C= 8.66, N= 7.74, 0= 8.57, Ne= 7.84, S= 6.87, 
Cl= 5.27, Ar= 6.23, Fe= 5.87. 

We do not attempt here to assign error bars to the abun- 
dances derived from the model, as this would require to per- 
fectly fit the models to various combinations of the extreme 
values of intensities allowed by the observations. This repre- 
sents a very time consuming task, because line intensities are 
not directly proportional to abundances. Fortunately, we do not 
need to go through such a lengthy procedure. What we wanted 
is to see whether we can obtain at least one photoionization 
model including all the classical ingredients which would fit 
all the observations, including the new infrared data, which put 
strong constraints on the temperature. The fact that we suc- 
ceeded argues against the presence of important temperature 
inhomogeneities in this object. 

5.3. Recombination lines and temperature 
fluctuation issues 

Let us now look at the recombination lines of C and O as 
predicted by our models. The C II A 4267 line is well repro- 
duced by our models, which indicates that the ionization cor- 
rection factors used by both Bresolin (2007) and Esteban et al. 
(2009) to derive the carbon abundance are in agreement with 
our model predictions. This, however, does not mean that the 
carbon abundance is correct, since one might expect a similar 
problem with the carbon recombination line as with the oxy- 
gen ones. Indeed, the sum of the intensities of multiplet 1 of 
O II (labelled Oil A 4651 in Table© is about 4 times smaller 
than the observed value (if we use Eq. 11 from Peimbert et 
al. 2005 to convert the intensity of O II A 4649 observed by 
Esteban et al. 2009 into the sum of the multiplet). This is be- 
cause the oxygen abundance in the models is fixed by fitting 
the observed intensities of the forbidden lines. Thus, the oxy- 
gen abundance discrepancy between the recombination lines 
and collisionally excited lines appears to be of a factor of 0.6 
That the models underpredict the value for O II A 465 1 is 
not surprising, since they are chemically homogeneous and do 
not experience sufficient temperature variations to explain the 
observed ++ /H + abundance discrepancy. 

We have made some experiments by enhancing the oxygen 
abundance so as to reproduce Oil A 465 1 . Oxygen cooling is 
then so important that it lowers the gas temperature to such a 
point where the [O III] /15007 line becomes severely underes- 
timated. In addition, if the oxygen abundance is enhanced by 
a factor of 3 of 4 with respect to the solar value, from what is 
known of stellar nucleosynthesis and galactic chemical evolu- 
tion, this should also be the case of other a-elements, such as 
neon, sulfur and argon. In this case, the intensitiy of the [S III] 
/U8.7/mi line becomes largely overestimated. 



5 This is larger than the value of 0.36 dex derived by Esteban et 
al. (2009) from their observations due to the conjunction of several 
causes such as the fact that our model attempts to fit the observation 
of Bresolin (2007) rather than those of Esteban et al. (2009), and that 
Esteban (private communication) did not employ the same procedure 
to estimate the sum of the multiplet. 



Our models also do not explain the Balmer jump tempera- 
ture of 5000 ± 800 K obtained by Bresolin (2007), which led to 
his t 2 = 0.06. However, it must be noted that the spectral res- 
olution of Bresolin 's observations was not sufficient to reach 
an accurate value of the level of the continuum close to the 
Balmer discontinuity as well as of the intensity of the high- 
order Balmer lines, so that the Balmer jump with respect to 
adjacent hydrogen lines could be easily overestimated by 20%, 
leading in that case to a Balmer jump temperature close to 7000 
K. The spectrum of Esteban et al. (2009) is an echelle spectrum 
with much higher spectral resolution but poor signal-to-noise 
ratio in the continuum, so unfortunately the Balmer jump could 
not be measured there. 

6. Conclusion 

In order to get more insight into the "temperature fluctuation" 
and the "abundance discrepancy" problems we have performed 
a detailed study of the giant H II region H 1013, for which ex- 
cellent spectroscopic data exist, and which has been previously 
diagnosed for the presence of temperature fluctuations. We 
have first completed the observational data base with our own 
mid-infrared measurements using the Spitzer Space Telescope, 
which give us important additional constraints, especially on 
the temperature structure. 

We then made a first estimate of the physical parameters 
(electron density and temperature) from different line ratios in- 
cluding the ones involving the infrared data using the PyNeb 
software (Luridiana et al. 2012). This indicated that tempera- 
ture fluctuations, even if present, would be much smaller than 
t 2 = 0.06 invoked by Bresolin (2007) from the Balmer jump 
temperature or t 2 = 0.04 from the ++ abundance discrepancy 
obtained by comparison of the recombination and collisionally 
excited lines from the high resolution spectrum of Esteban et 
al. (2009). 

For a deeper insight, one needs to carry out a detailed pho- 
toionization analysis in which all the processes that affect the 
line intensities as well as the observing apertures are taken into 
account. We performed such a study, using Starburst99 and the 
Cloudy _3D package, trying to find a photoionization model of 
H 1013 that is able to reproduce all the observational data (sur- 
face brightness distribution, line ratios and stellar features). 

We were able to produce an almost satisfactory model, as- 
suming a volume filling factor of 0.005. Since the use of a fill- 
ing factor is just a zero-order approach, we checked the effect 
of adding a low-density medium to the modelled clumps or fil- 
aments. This resulted in an even better fit to the observations. 
Of course, as in any study which requires the inclusion of a vol- 
ume filling factor or - more realistically - with a large volume 
of diffuse gas bathing the filaments, we are left with the ques- 
tion of what continuously generates those filaments, since their 
lifetime must be short, as they are not in pressure equilibrium 
with the surrounding diffuse matter. 

The fact that we are able to produce photoionization mod- 
els that perfectly fit most observational constraints - which are 
strong and numerous - argues against the presence of signifi- 
cant temperature fluctuations in H 1013, contrary to what has 
been advocated previously for this object. The oxygen abun- 
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dance of our best model is 12 + log O/H = 8.57, as opposed to 
the values of 8.73 and 8.93 advocated by Esteban et al. (2009) 
and Bresolin (2007), respectively. 

However, our model is not able to reproduce the intensi- 
ties of the oxygen recombination lines, which is one of the 
big things that remain to be explained. Our model also does 
not reproduce the very low B aimer jump temperature inferred 
by Bresolin (2007) from his data, but we have argued that his 
value might be in error, possibly due to insufficient spectral 
resolution of the observations for such a purpose. A new obser- 
vational study, focussing on the B aimer jump temperature in 
this nebula, would obviously be very important for our case. 

This paper can by no means be considered a solution to 
the nagging "temperature fluctuation" and "abundance discrep- 
ancy" problems, since it studies just one object and does not 
even provide a complete understanding of the observational 
data. But it serves as a warning against too general statements 
about the physical conditions in H II regions and too hasty con- 
clusions drawn from their spectra. A promising solution is of 
course the deviation from a Maxwellian distribution of electron 
energies, proposed by Nicholls et al. (2012). Work in this direc- 
tion is needed to show quantitatively that such circumstances i) 
can occur in H II regions and ii) are able to produce photoion- 
ization models that reproduce a wide range of observational 
constraints, as we have attempted. 
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